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We propose that the properties of glass transition can be understood on the basis of elastic 
waves. Elastic waves originating from atomic jumps in a liquid propagate local expansion due 
to the anharmonicity of interatomic potential. This creates dynamic compressive stress, which 
increases the activation barrier for other events in a liquid. The non-trivial point is that the range 
of propagation of high-frequency elastic waves, d e \, increases with liquid relaxation time r. A self- 
consistent calculation shows that this increase gives the Vogel-Fulcher-Tammann (VFT) law. In 
the proposed theory, we discuss the origin of two dynamic crossovers in a liquid: 1) the crossover 
from exponential to non-exponential and from Arrhenius to VFT relaxation at high temperature 
and 2) the crossover from the VFT to a more Arrhenius-like relaxation at low temperature. The 
corresponding values of r at the two crossovers are in quantitative parameter-free agreement with 
experiments. The origin of the second crossover allows us to reconcile the ongoing controversy 
surrounding the possible divergence of r. The crossover to Arrhenius relaxation universally takes 
place when d e i reaches system size, thus avoiding divergence and associated theoretical complications 
such as identifying the nature of the phase transition and the second phase itself. Finally, we discuss 
the effect of volume on r and the origin of liquid fragility. 



INTRODUCTION 

The problem of glass transition has been widely dis- 
cussed [l], LH, S 0, S IE 0> an d nas been considered as 
one of the deepest and most interesting challenges in 
physics 0. As widely perceived, a glass transition theory 
should provide a consistent explanation of several uni- 
versal properties of liquids which set in on lowering the 
temperature, including the physical origin of the Vogel- 
Fulcher-Tammann (VFT) law, slow non-exponential re- 
laxation and dynamic crossovers [J, LH, LH, EL IE IE 0, S 
S LH LH, LH, Il4j . The most widely studied property is 
the unusual behaviour of liquid relaxation time, r. On 
lowering the temperature, t is almost never Arrhenius, 

but follows the VFT law: t = Tq exp ^ T ^ To ) > where A 

and T are constants 0, LI, i, LI i, @, ■ 

As recently reviewed the quest to understand the 
origin of the VFT law and other anomalous features of 
glass transition has resulted in the development of many 
theories and models, which discuss different parameters 
that control glass transition: free volume, entropy, energy 
landscape, mode coupling and others. However, there 
is no agreement as to what physical parameter governs 
glass transition [2|. For this and other reasons, it has 
been proposed that glass transition remains a mystery, 
with no simple picture emerging [l|, LH, 0, HI LH, LH, 0] • 

Existing theories are often elaborate, and approach 
glass transition as an outstanding phenomenon that re- 
quires novel or special ideas and mechanisms Q. Yet we 
feel that one should be able to describe the process of 
cooling a liquid to a glass using familiar physical con- 
cepts which may, however, operate in a non-trivial and 
unexpected way. We suggest that a property relevant to 
glass transition is elasticity because a glass differs from 
a liquid only by its ability to support static shear stress. 



Hence, we approach the problem by asking whether glass 
transition can be understood on the basis of liquid elastic 
properties. 

Elastic approaches to glass transition were discussed 
previously (see, e.g., Ref. @, 



15|, ll6|), but the prob- 



lem of explaining glass transition from the first principles 
remains. Consequently, there is no microscopic under- 
standing of the origin of the VFT law, slow relaxation, 
dynamic crossovers and other effects of glass transition. 

Recently, we proposed that several important proper- 
ties of glass transition can be understood on the basis of 
elastic waves in a liquid 17 , 1|| . We considered the case 
when a liquid is perturbed and relaxes to equilibrium. 
However, it is important to consider elastic waves that 
originate from local atomic jumps in an equilibrium liq- 
uid. The main question is why and how these waves can 
result in the slowing down of liquid dynamics, the VFT 
law and dynamic crossovers. 

In this paper, we develop and extend our approach to 
glass transition. Considering several anomalous, yet uni- 
versal relaxation laws discussed above, we explore non- 
trivial and unusual ways in which liquid elastic prop- 
erties may emerge during glass transition. We analyze 
elastic waves originating from atomic jumps in a liquid, 
and find that their effect is to create a dynamic compres- 
sive stress which slows down relaxation of other events. 
The non-trivial point is that the range of propagation of 
these waves increases with liquid relaxation time. A self- 
consistent calculation shows that this increase gives the 
VFT law. In the proposed theory, we discuss the origin 
of two dynamic crossovers in a liquid, the absence of di- 
vergence of r at To, the effect of volume on r and the 
origin of liquid fragility. 
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ELASTIC INTERACTION BETWEEN LOCAL 
RELAXATION EVENTS 

Unlike in solids, atoms in liquids are not fixed in space, 
but are constantly rearranging. This gives liquid flow. 
Each flow event is a jump of an atom from its surround- 
ing "cage" , accompanied by large-scale rearrangement of 
the cage atoms. We call this process a local relaxation 
event (LRE). A LRE lasts on the order of Debye vibra- 
tion period r wO.l ps. 

There are two known basic properties of LREs. The 
first property concerns liquid relaxation time, r. r was 
phenomenologically introduced by Maxwell in the vis- 
coelastic picture of flow as t = ri/Goo, where rj is liq- 
uid viscosity and Goo is the instantaneous shear modulus 
[lit . Frenkel offered microscopic interpretation of r as 
the time between LREs at one point in space in a liquid 
[20l | . At high temperature, r w To. When r increases to 
t w 10 3 s at glass transition temperature T g , a liquid, by 
convention, forms a glass 0] ■ The second property is that 
a LRE requires an increase of local volume. As widely 
discussed 0, [2(| , the activation barrier for a LRE at con- 
stant cage volume is very large due to strong short-range 
interatomic repulsions. Hence, atoms in the cage need 
to increase its volume in order to allow for the escape of 
the central atom (see Fig la). In doing so, work is per- 
formed to deform the surrounding liquid. This probes liq- 
uid elasticity. The work against the elastic force is equal 
to the activation barrier for a LRE, U 0, H3]- This bar- 
rier is surmounted by temperature fluctuations, so that 
t = T exp(C//T) (fc B = l)H- 

An important insight into glass transition comes from 
the realization that LREs interact elastically, as we have 
recently proposed [13, 03 • A LRE involves restructuring 
of the cage that involves large-amplitude atomic motions 
of about of 1-2 A (see Fig. la). On the very short time 
scale of a LRE (when To < t), the surrounding liquid 
can be viewed as an elastic medium (2(|. Therefore, the 
large atomic motion from a LRE elastically deforms the 
surrounding liquid, inducing elastic waves. Because their 
wavelengths are on the order of interatomic separations, 
the frequency of these waves, oj, is on the order of De- 
bye frequency. This means that in almost entire range 
of t that is relevant for glass transition, u> > 1/t holds 
true. As discussed by Frenkel, high-frequency ui > 1/t 
waves are propagating in a liquid as in a solid [20] . The 
existence of propagating high-frequency waves in liquids 
is now firmly established: the conclusion from numer- 
ous experiments is that liquids support vibrational modes 
which extend down to wavelengths comparable to inter- 



21] 



atomic separations, similar to phonons in solids 

The propagating high-frequency waves from a LRE de- 
form the cages around other atoms in a liquid. This af- 
fects their relaxation because, as discussed above, the 
jump of an atom depends on its cage. Therefore, LREs 
interact via the elastic waves they induce. Elastic interac- 
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FIG. 1: (a) Large-scale cage restructuring due to the atomic 
jump induces a propagating high-frequency wave. This wave 
propagates volume expansion due to anharmonicity; (b) As a 
result of the arrival of compressive wave fronts, atoms in the 
central cage are under dynamic compressive stress. 



tion between LREs is the physical origin of cooperativity 
of relaxation in a liquid, whose physical origin has been 
much discussed, but remained unclear EL The 
key issue is the range of this interaction. 

Lets consider how a solid-like elastic wave is affected by 
LREs in a liquid. Because t sets the period of structural 
rearrangements in a liquid, it defines the time of decay 
of induced static shear stress [20|, 122J as well as of high- 
frequency (uj > 1/t) solid- like propagating waves [23(. 
Note that t is similar for high-frequency shear and lon- 
gitudinal stress [23( 1 . If c is the speed of sound, d e \ = ct 
gives the length of stress decay: 



d c i = ct (1) 

At the microscopic level, d e \ = ct originates as follows. 
A high-frequency u) > 1/t wave induced by a LRE prop- 
agates as in a solid until a remote LRE takes place at 
the wave front. At this point, atoms in the wave front do 
not pass the oscillations further on as in a solid with con- 
stant structure, but are involved in large scale structural 
fluctuation due to the LRE motion. Hence, the solid-like 
elastic wave is absorbed by the motion of the liquid-like 
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atomic cluster due to the LRE. Suppose the remote LRE 
meets the wave front distance d e \ away from the original 
LRE. d e i is defined from the equality of the wave travel 
time, d e \/c, and the time at which the remote LRE takes 
place at point d e \. The latter time is given by r because 
microscopically, r is the average time between LREs at 
one point in space, and we find d c \ = ct as before. 

Hence, d Q \ = ct is the distance over which a high- 
frequency wave propagates in a liquid as it would in a 
solid without its structure being modified by LREs. In 
an ideal crystal with infinite r, Eq. |lj gives the infinite 
range of wave propagation, as expected. 

It is interesting to note that d c \ = ct is in agreement 
with Frenkel's theory of viscoclastic relaxation in a liquid. 
In this theory, waves in a liquid decay with distance x as 
oc exp(— x/d). Using the condition lot > 1 explicitly, 
Frenkel's theory gives d » ct. This is discussed in detail 
in the Appendix. 

We emphasize that similar to Frenkel's theory, our 
derivation of d e \ = ct is also based on the condition 
lot > 1, albeit implicitly. Indeed, we approached a liq- 
uid from the solid (elastic) phase, and considered how 
solid-like waves are affected by LREs in a liquid. This 
approach is therefore based on the assertion that solid- 
like elastic equilibrium exists in a liquid, which is the case 
for high-frequency to > 1/r waves. 

The non-trivial point here is that d c \ = ct increases 
with r. This is directly opposite to the usual decay of 
hydrodynamic waves, whose propagation range varies as 
1/r (see also, Appendix). Crucially, this is because the 
considered solid-like elastic regime of wave propagation 
(lot > 1) is markedly different from the commonly dis- 
cussed hydrodynamic regime (lot < 1). In the latter, 
LREs are frequent enough to eliminate the state of elas- 
tic equilibrium, and establish hydrodynamic equilibrium 
instead. The resulting equations of motion are those of 
hydrodynamics and viscous flow 24j . Interestingly, these 
equations are widely used to describe liquid dynamics 
and glass transition, yet they are not applicable to our 
approach to glass transition based on high-frequency in- 
teractions. The presence of these interactions makes our 
approach essentially non-hydrodynamic, but elastic in- 
stead. 

d e i can be called liquid elasticity length, because it de- 
fines the range over which two LREs interact with each 
other via induced high-frequency elastic waves. Impor- 
tantly, d c \ = ct increases on lowering the temperature 
because r increases. We propose that this is the key to 
the problem of glass transition. 

We finish this section with a comment regarding the 
generality of our discussion. Depending on liquid struc- 
ture and interactions, LREs may take different form. For 
example, a LRE in covalent network liquids may involve 
bond switching from under- to over-coordinated states, 
whereas in spherically symmetric systems (e. g., ionic or 
metallic) it may resemble an illustration in Fig. 1. Hence, 



the way in which high-frequency waves are generated by 
LREs may be system-specific. However, the interaction 
of LREs via induced elastic waves is general, and should 
apply to all liquids undergoing glass transition, including 
covalent, ionic, molecular, metallic, polymeric and oth- 
ers. Consequently, we expect that the VFT law, dynamic 
crossovers and other effects of glass transition discussed 
below can be understood on the basis of elastic interac- 
tions and d e i. 



DYNAMIC COMPRESSIVE STRESS 

Lets consider the nature of LRE-induced waves in more 
detail. In between LREs, atoms in a given local region vi- 
brate with small amplitudes as in a solid (glass) [20(, and 
a harmonic approximation can be applied to the same ex- 
tent as to the solid phase. As discussed above, a LRE 
involves large-scale atomic motions of 1-2 A. This consid- 
erably widens the distribution of interatomic separations 
Ar. At large Ar, harmonic approximation no longer ap- 
plies, and the potential anharmonicity becomes impor- 
tant. In particular, the decrease of Ar results in short- 
range repulsion that is always stronger than the attrac- 
tion due to the same increase of Ar. Because expansion 
carries smaller energy penalty, large-amplitude motion of 
atoms involved in a LRE results in short-lived expansion 
of local volume around the LRE. In a simplified picture, 
this process can be thought of as the appearance of "hot" 
local regions in a liquid due to LRE motion and associ- 
ated local thermal expansion. 

Local volume expansion is propagated away from a 
LRE by high-frequency elastic waves discussed above. 
We note that when a sphere expands in a static elas- 
tic medium, no compression takes place at any point. 
Instead, the system expands by the amount equal to the 
increase of the sphere volume resulting in a pure 
shear deformation. The strain components u from an 
expanding sphere (noting that u — >0 as r — > oo) are 
u rr = —2b/r 3 , ugg = it00 = b/r 3 [22|, giving pure shear 
uu = 0. As a result, the energy to statically expand 
the sphere of radius r by amount Ar depends on shear 
modulus G only: E = 8nGrAr 2 0. 

Unlike in the static case, there is compression at the 
front of the expanding wave. At time t, the wave front 
causes an outward displacement of atoms on the sphere 
of radius ct. This takes place during time To, i.e. very 
fast, because the wave is due to the remote LRE that lasts 
approximately To. Until this displacement causes the mo- 
tion of the next concentric sphere with radius of ct + x, 
where x is on the order of the interatomic separation, 
there exists a brief compressive stress in a layer of thick- 
ness x. This stress exists over time approximately equal 
to the period of vibrations of the system of two atoms, 
or To, i.e. is very short. As the wave propagates further, 
dynamic compression between the two layers disappears, 
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and shear deformation is established as discussed above. 

Lets now fix an atomic cage in the centre in Fig. lb. 
As the front of the wave propagating volume expansion 
from a remote LRE arrives at the centre, it causes a brief 
compressive stress between the cage boundary atoms and 
the atoms in the next layer that include the central atom 
(see Fig. lb). Hence, it puts the cage atoms under dy- 
namic compressive stress (DCS). Because, as discussed 
above, the jump of the central atom requires cage expan- 
sion, the compressive wave fronts arriving at the centre 
result in more work needed to expand the cage. This 
increases the activation barrier for the central LRE, U, 
and slows down liquid dynamics. 

We emphasize that DCS in a liquid is created by new 
elastic waves, which are notably absent in solids, and 
which propagate local volume expansions from the hot 
(in the sense discussed above) local regions. This process 
can be compared to pressure waves generated by local 
laser heating of a solid. In this case, the energy that 
heats up local regions and creates these pressure waves is 
external. In a liquid, the energy to create DCS is internal, 
i.e. it is the liquid's thermal energy. Here, as temperature 
increases and the liquid state is approached from the solid 
phase, LREs appear as a new type of "hot" local motion. 
This motion gives rise to a new set of expanding elastic 
waves, and DCS emerges as a result. 



DERIVATION OF THE VFT LAW 

The stage is now set for the calculation of U. U is 
equal to the total work required to expand the cag e by 
the amount required for a LRE to take place |20j. U 
can be written as U = Uq + U\. Here, Uq is the high- 
temperature activation barrier that depends on liquid in- 
teratomic forces and structure, but not on the elastic in- 
teractions with other LREs, i.e. is non-cooperative, or in- 
trinsic. U\ represents the interaction (cooperative) term. 
U\ is equal to the additional work to expand the cage due 
to the arriving compressive wave fronts from other LREs, 
and is set dynamically. Let q be the increase of the cage 
volume required for a LRE to take place, q « a 3 , where a 
is the interatomic separation of about 1 A. As discussed 
above, DCS exists in the arriving wave front, in a layer 
of thickness a. Lets consider DCS created by a remote 
LRE i distance r away from the centre in Figure lb. If 
Pi(r) is the value of DCS at the centre, the remote LRE 

JVt 

contributes qpi(r) to U\. Then, U\ — ^ qPi(r), where 

i=i 

N T is the number of compressive wave fronts that pass 
through the centre during time t, and U reads: 



U 



i=l 



(2) 



According to Eq. (JTj), the sum in Eq. ([2J includes 
elastic waves from LREs inside the sphere of radius d e \ — 
ct. Importantly, each of the local relaxing regions inside 
this sphere contributes to the sum once. This is due 
to two reasons. First, because the central event relaxes 
during time r, waves from all the events located distance 
ct away from the centre have enough time to propagate 
to the central point. Second, because a remote event 
also relaxes during time t, it contributes only one wave 
during the time of relaxation of the central event. Hence, 
the introduced length d c \ = ct is self-consistent in that 
it accounts for the dynamical nature of LREs as well as 
for the wave dissipation. 

Therefore, Eq. ([2|) can be written as 



U = U + 



<! I 4irpr pi(r)dr 

d /2 



(3) 



where p = is the density of local relaxing regions and 

do is the region diameter of about 10 A. 

Pi(r) decreases with r. Recall that the short-lived DCS 
from a LRE originates at time t in a thin layer of thick- 
ness x, when the outward motion of atoms takes place in 
a sphere of radius ct, but not yet in a sphere of radius 
ct + x. The strain and stress on the sphere of radius ct 
decay as in elastic medium. This is because the waves 
under consideration are of high frequency (u> > 1 /r) and 
therefore propagate in elastic equilibrium as discussed 
above. In elastic medium, the strain on an expanded 
sphere (i.e. on the sphere of radius ct) decays asuoc^ 
|22j |. This strain creates pi which, therefore, also decays 
as Pi(r) oc 4f- Let po be the value of pi at the cage 
boundary, distance do/2 away from the centre of a relax- 
ing region. Then, Pi(r) — p (|^) , and Eq. becomes 



U = Un 



3qpo hi 



24] 
da 



(4) 



According to Eq. (g]), U increases with d c \. Because 
d e \ = ct itself increases with r and hence, with U, U is 
defined self-consistently. Combining d e \ = ct with r = 
Tq exp(U/T) and noting that a « ctq, we write 



d e i = a exp I — 
Putting this in Eq. (g]), we find 



(5) 



U 



AT 
T-T n 



(6) 



where T = 3qp and A = U - T In ^. From Eq. 
the VFT law follows. 
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In this theory, the origin of the VFT law is the increase 
of d e i on lowering the temperature (see Eqs. (|TJ) and 
Q). This increase results in a larger number of LREs 
that elastically interact with a given event, increasing its 
activation barrier. The transition from the VFT to the 
Arrhcnius law takes place in the limit of small d e \ at high 
temperature. In this case, Eq. (|4]) gives U = Uq, i.e. U 
becomes non-cooperative and temperature-independent. 
This gives Arrhcnius relaxation: r = tq exp([/n/T). 

To = 3qpo can be roughly estimated by recalling that 
Po is created by local volume increase due to anharmonic- 
ity. Because this increase gives rise to macroscopic ther- 
mal expansion, po can be estimated as po = aBT m , where 
B is bulk modulus, a is the coefficient of thermal expan- 
sion and T m is the temperature of local motion during a 
LRE, which is on the order of melting temperature. Tak- 
ing q ps 1 A 3 and typical liquid values of B ps 10 GPa 
and a = 1CT 4 - 1CT 3 K" 1 , we find T = (0.1 - l)T m , in 
ordcr-of-magnitude agreement with experimental values 
(§]. Hence, To in the derived VFT law is related to liquid 
parameters that are physically sensible. 



EFFECT OF VOLUME 

In addition to temperature, our theory also predicts 
the dependence of t on volume or pressure. According to 
recent experiments, r is a function of both temperature 
and volume, although temperature has a larger overall 
effect [25[. In our theory, pressure has two effects: first, 



q becomes larger by the amount equal to the decrease of 
the cage volume due to external pressure P, so that the 
new q' = q + g c |, where q c is the initial cage volume. 
Second, the activation barrier U in Eq. and hence 
the VFT parameter A, acquire an additional term q'P 
due to the increased work against the external pressure. 

An interesting prediction is the dependence of To on 
pressure: because To = 3qpo (see Eq. (JBJ), its value 
under pressure, Tq, is Tq = 3poq' = 3po (q + q c ^) = 
T + 3p Q q c ^. For small P, B increases linearly with pres- 
sure: B = Bq + CP, where B is zero-pressure bulk mod- 
ulus and C is a constant, giving Tq = T + 3poq c 



ic Bo+CP- 

Then, for small Tq = T + Zp Q q c ^ (l - C^) • This 
behaviour is observed in the experiments: at small P, T 
increases linearly with P. with a negative quadratic term 
appearing at higher P 26l. 21\. 



DYNAMIC CROSSOVERS 

We now discuss how our theory explains the origin 
of dynamic crossovers, the important problem of glass 
transition 0, 0, S i, 0, S M, HH, 0, 111 H]. Exper- 
imentally, there are two dynamic crossovers in a liquid. 
The first crossover is at high temperature, and marks the 



transition from exponential to slow stretched-exponential 
dynamics and from Arrhenius to VFT relaxation. The 
physical origin of this crossover has remained one of the 
central open questions in the area of glass transition 
0, B B 0, H @LFor various liquids, r at the crossover is 

t « i -30 ps [13, El- 

In our theory, the origin of this crossover is understood 
as follows. At high temperature when t — To, Eq. (fT]) 
gives G? G i = cto = a ~ 1 A. Therefore, at high tempera- 
ture elastic waves from LREs do not propagate beyond 
the nearest-neighbour distance. This means that d e \ ps 1 
A is shorter than the distance between two neighbouring 
LREs, or two adjacent molecular cages, d m , which varies 
from d m ps 10 A in small- molecule to d m ps 100 A in 
large-molecule liquids. Because d e \ < d m , LREs do not 
elastically interact and, therefore, relax as independent, 
resulting in exponential and Arrhcnius relaxation. On 
the other hand, when d e \ increases to d m on lowering the 
temperature, LREs are no longer independent, but start 
interacting via the induced elastic waves. This interac- 
tion gives the VFT law as discussed above and, as we 
recently showed [17| , stretched-exponential relaxation. 

Therefore, the first crossover corresponds to d c \ = d m . 
From Eq. (UJ, r at the first crossover, t\, is 



Tl 



= T , 



(7) 



giving Ti Ri (10 — 100)to ps 1 — 10 ps, consistent with the 
experimental results. 

The second dynamic crossover is at low temperature, 
and marks another qualitative change in liquid dynamics 
11, 13, 13, | . The important change is the crossover 
from the VFT law to a more Arrhenius behaviour. Start- 
ing from the early work in Ref . [28| , it was realized that 
at low temperature the VFT law predicts viscosity and 
U that are larger than those experimentally measured. 
Moreover, the experimental U at low temperature be- 
comes temperature-independent, contrary to its contin- 
uous increase predicted by the VFT law. The crossover 
from the VFT law to a more Arrhenius form at low tem- 
perature has now been established in a large number 
of glass- forming liquids 11 , [l3|, EH • The origin of this 



crossover is not understood at present. 

In our theory, the origin of this crossover is as fol- 
lows. When d c \ — L, where L is system size, all LREs in 
the system elastically interact. This gives temperature- 
independent U oc ln(T) in Eq. Q. Hence, U can not 
increase due to the increase of d c \ on lowering the tem- 
perature, but due to other effects only (e.g., density in- 
crease). As a result, t crosses over to a more Arrhenius 
form. 

Hence, the second crossover corresponds to d c \ = L. 
From Eq. ([TJ) , r at the second crossover, T2, is 
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L L 

T 2 = — = —T 

c a 



(8) 



Taking a typical experimental value of L in the range of 
0.1-10 mm, we find r 2 w 10~ 7 — 10~ 5 s. This is in good 
agreement with r of the crossover from the VFT to a 
more Arrhenius relaxation seen experimentally [14| . In 
addition, r 2 agrees well with the experimental r at which 
other important liquid properties show a crossover and 
undergo qualitative changes [12l | . 

Therefore, the two dynamic crossovers originate in 
our theory in a simple and physically transparent way. 
Derived solely from the definition of d e \ in Eq. {T|), 
ti = d m /c and r 2 — L/c give good agreement with the 
experiments, without using adjustable parameters. This 
lends support to the proposed theory of glass transition. 

Importantly, the origin of the second dynamic 
crossover in this theory allows us to reconcile what is 
perhaps the main ongoing controversy surrounding glass 
transition, that of divergence [l], 0, 0, HI, @] ■ Formally, r 
in the VFT law diverges at To. Although To is always 
smaller than T g , it is natural to ask what the physical 
significance of To is. Starting from early theories of glass 
transition, To has been associated with a phase transi- 
tion into a state of zero configurational entropy termed 
"ideal glass" @, 0, 0, @]. Confounded by a number of 
problems (H, this approach was followed by other diver- 
gence scenarios of glass transition that were based on 
the existence of an underlying or avoided phase transi- 
tion [B[ . However, the nature of the phase transition and 
the second phase itself remain unclear, primarily because 
no long-range order or any other appreciable structural 
changes appear on cooling 0, @|. Moreover, as recently 
emphasized, there are no experimentally observed signs 
of divergence [B| . These issues continue to fuel the current 
debate about whether divergence exists and if so, what it 
means. The problem is often formulated as whether glass 
transition is a thermodynamic or dynamic phenomenon 

@,aaa. 

In our theory, this controversy is reconciled as follows. 
Combining Eqs. ([5]) and Eq. ©, we write 



del = a exp 



.4 



T-T 



(9) 



When T approaches To, d c \ quickly exceeds any finite 
size of the system. At this point, r crosses over to more 
Arrhenius as discussed above, moving the divergence to 
zero temperature. Therefore, r does not diverge at To. 

We note that in addition to a seeming divergence of r 
at To, approaches to glass transition based on thermody- 
namics and phase transitions have also been stimulated 
by the experimental changes of heat capacity, compress- 
ibility and other properties at T g 0, 0, S| ■ As we recently 
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FIG. 2: r as a function of T g /T for Si02 (°) and o-terphenyl 
(o). t is calculated from r\ = GooT, where Goo ~ 10 GPa. 
The data of r\ are from Ref. [7|. 



discussed [29[ , these changes can be understood as a nat- 
ural signature of T g insofar as T g is defined by the freezing 
of LREs at the experimental time scale. In other words, 
the observed changes are related to a liquid falling out of 
equilibrium at T g , rather than to thermodynamics or a 
phase transition [291 ] . 



FRAGILITY 

Our theory readily explains the origin of liquid fragility 
[E 0,0, 00 , a widely debated subject in the area of glass 
transition. Such an explanation was offered in our recent 
paper [18j |. which we briefly recall below. 

Lets consider two extreme cases of "strong" Si0 2 and 
"fragile" o-terphenyl (see Fig. (2)) and calculate d c \ at 
the highest measured temperature, T/j. From Fig. (2), 
r(Th) is approximately 10~ 7 s and 10 -12 s for Si0 2 and o- 
tcrphcnyl, respectively. Taking c w 1000 m/s, d e \(Th) is 
about 0.1 mm for Si0 2 and 1 nm for o-terphenyl. Hence, 
del for o-tcrphenyl grows from microscopic to macro- 
scopic values on lowering the temperature. This gives 
large increase of cooperativity and U in Eq. (f5|), i.e. 
fragile behaviour. On the other hand, d e \ for Si0 2 ap- 
proaches system size at Ti already, leaving little room 
for the cooperativity and U to grow on lowering the tem- 
perature (see Eq. ([!])). This gives strong behaviour. 

In this theory, liquids that have large intrinsic activa- 
tion barrier and large viscosity (r) also have large d e \, 
even at high temperature. This gives little room for U to 
grow, resulting in strong behaviour. On the other hand, 
liquids that have small viscosity at high temperature and, 
consequently, small d c \, show large increase of U on low- 
ering the temperature, i.e. are fragile. An interesting 
prediction of our theory is that if the measurements are 
extended to higher temperature so that r and d e i are 
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small, liquids will become more fragile. Note that there 
is no data for strong liquids in the range of small r 
(see also, Fig. (2)) due to high melting points, and the 
evidence for Arrhenius behaviour comes only from the 
range where r is large. 

In our theory, fragility can be quantitatively related 
to other system properties. Fragility is quantified by pa- 
rameter D = 7] ; the larger D the smaller fragility. In 

our theory, A = U - T In |a (see Eq. ^ and below), 
giving D = ^ — In 7^. Therefore, fragility is predicted 
to decrease with the scaled high-temperature activation 
barrier Uq. This is consistent with the experimental data 
showing the decrease of fragility with ^ [lit . 

COMMENT ON DCS 

We make two remarks about DCS. First, if the fre- 
quency of compressive wave fronts arriving at the centre 
in Fig. lb is large enough, the atoms in the central cage 
may, at first glance, appear to be under an effective static 
stress. The necessary condition for this is that the aver- 
age time difference between any two compressive arriving 
wave fronts, At, is shorter than the minimal time set by 
the elementary vibration period tq: At < tq. As dis- 
cussed above, the central region is affected by all events 
inside the sphere of radius d c \ = ct, each contributing 
once. Hence, there are N T — r^T^p contributing lo- 
cal relaxing regions inside the sphere of radius cr. The 
central event relaxes during time t, hence the average 
time between the wave fronts arriving at the centre is 

(d /2)^ 

At = -jq- = 1 °J T J . Because c ~ a/ To, the condition 

3 

At < t gives r > t (1^-) 2 . Denoting 

Tram = T (J^J ( 10 ) 

we find that r m i„ ~1 ps, i.e. is very short. 

As discussed in the previous section, d c \ exceeds the 
experimental system size L above T g . In this case, we 
should substitute ct for L, and the condition At < tq 

gives r < r ( jfc j . Denoting 




and using a typical value of L of 1 mm, we find T max ~ 
10 5 s. Because r max » t(T 9 ), the condition At < t is 
also satisfied when d c \ > L for a macroscopic system. 

Therefore, the necessary condition for the cage atoms 
to be under the static stress (At < To) is fulfilled in a very 
wide range of r, including glass transformation range. 
However, the actual value of this stress, p s , vanishes for 



t > To. Indeed, the contributions from arriving compres- 
sive wave fronts to the static stress can be summed on the 
time scale of To only, for two reasons. First, wave fronts 
from remote LREs add up only if they arrive simultane- 
ously, or in practice during the elementary time period 
To. Second, the duration of the passing wave front is also 
To, because the lifetime of the remote LRE that creates 
the wave is about to. Therefore, p s can be calculated as 

Ps = Pt J 4:Trr 2 pi(r)dr (12) 

d /2 

where p To is the density of those remote local regions that 
give rise to the waves passing through the centre during 
time To- 

The number of the wave fronts that pass through the 
centre during time To is iV ro = ft. Combining this with 
At = from above gives N T() = ^N T . Hence, p T() — 
^p, where p — -|j is the density of local relaxing regions 
introduced in Eq. ([3]). Using p TQ in Eq. p2[) and recalling 
that d i — ct = and pi(r) = p (|a ) 3 , we write 

p.=3po-fln-+ln^ > ) (13) 
t V T o do J 

Therefore, p s vanishes for r 3> tq, i.e. for r at which 
del and elastic interactions between LREs become appre- 
ciable in the first place (see the previous section). As a 
result, DCS does not contribute to liquid internal pres- 
sure or equilibrium volume. Physically, this is because 
the volume of local regions that contribute to p s is sparse. 

We note that the vanishing of p s is in contrast to the 
behaviour of U which increases with r and d e \ (see Eq. 
dU)). U is given by the total work to be performed against 
cage expansion [20| ■ Because each LRE compresses the 
central cage, this work, being an extensive quantity, is 
the sum over all pulses that arrive during time r from 
within the sphere of radius d c \. Therefore, U increases 
with the number of contributing LREs and d e \. 

The second remark about DCS concerns Eq. |(5J). Af- 
ter the initial expansion due to a LRE, a remote local 
region relaxes back to its original equilibrium volume, 
sending a dilatational wave to the centre. However, the 
dynamic dilatational stress, pa, and the dynamic com- 
pressive stress, p, do not oscillate around an equilibrium 
value as in a harmonic wave, but are set by processes that 
operate on different time scales (or, equivalently, length 
scales, see below), making propagating density variations 
essentially anharmonic. Recall that p is due to fast ex- 
pansion originating from anharmonic forces that initially 
appear between several (3-5) "hot" atoms of the remote 
LRE during the elementary vibration period t ~ 0.1 
ps. As the expanding wave propagates, fast dynamic 
compression is created between the two adjacent atomic 
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layers when one layer is already displaced during tq but 
the other is not. 

On the other hand, dilatation that sets pd and origi- 
nates from cooling and contraction of the remote relaxing 
region is a slower process for the following reason. Let 
td be the time during which pd exists, td is set by the 
time during which increased interatomic separations in 
the expanded remote region, Ar, return to their original 
(pre-LRE) values, Ar . This process involves the number 
of atoms that is necessarily larger than the number of ini- 
tially hot atoms, because these hot atoms interact with 
their neighbours. As a result of this interaction, large 
Ar between the initially hot atoms are distributed over 
a larger region. This is accompanied by a gradual reduc- 
tion of Ar. As thermalization of the local region proceeds 
and Ar = Aro is established for all atoms involved, the 
dilatation of the relaxing region is complete. Therefore, 
td is set by the time of thermalization of atoms affected 
by the LRE motion. This thermalization involves at least 
the nearest neighbours of the initially hot atoms, hence 
td is the time of thermalization of the region whose size 
is at least equal to the cage size do- Therefore, the lower 
limit of td can be estimated as « 1 ps, because ther- 
malization can not proceed faster than the phonon speed, 
and we find td >1 ps. Hence, td^$> t . 

Interestingly, a similar effect is observed in molecular 
dynamics simulations of radiation damage. Here, a hot 
radiation cascade created by an energetic ion elastically 
deforms the surrounding lattice. It is found that fast 
initial expansion of the lattice is followed by considerably 
slower contraction due to finite thermal conductivity [301 ] . 

td ^ tq means that pd <C P- A mechanical analogy 
of this effect is compression of a spring due to a fast 
compressive force applied to one end, followed by slow 
motion in the opposite direction during which the spring 
length hardly changes. 

— can be related to t 1 as follows. Lets consider that 

p t d 

the boundary atom in the central cage is pushed towards 
the central atom by force / due to compressive stress p, 
followed by the motion in the opposite direction caused 
by force fd due to dilatational stress pd that arrives later. 
Let I be forward and reverse displacement of the bound- 
ary atom due to p and pd, respectively. Assuming con- 
stant acceleration, / = and fd = ^%r, where m is the 

d 

2 



atomic mass. Then, = J jT = ff^J ■ Hence, pd <C p 
because tq -C t<j. 

— can also be written as the ratio of the short-range 

and medium-range order distances, j^. Using td > ™ 

c = f = (^) 2 < ( a.) 2 . Hence, p d « p be- 

cause a <C do- Physically, this result has the following 
meaning. A statement equivalent to tq <C t d is the asser- 
tion that the wavelength of the initial compressive wave 
is comparable with interatomic separations in the central 
cage (because this wave is created by fast expansion of 



the remote cage), whereas the wavelength of the follow- 
ing dilatational wave is larger than the cage size due to 
slow contraction. Hence, interatomic separations in the 
central cage decrease in the first case, whereas the cage 
moves as a whole in the second, resulting in pd <C p. 

From pd <C p, qpd <C qp follows, i.e. the negative work 
due to pd is small, and can be ignored in Eq. I[2|). 

SUMMARY 

In summary, we proposed that the properties of glass 
transition can be understood on the basis of elastic waves 
in a liquid. Elastic waves originating from atomic jumps 
in a liquid create a dynamic compressive stress, which 
slows down relaxation. The increase of d e \ on lowering 
the temperature gives the VFT law. In addition to tem- 
perature, we also predicted the effect of volume on r. In 
the proposed theory, we discussed the origin of dynamic 
crossovers, the absence of divergence of r at To, the effect 
of volume on t and the origin of liquid fragility. 

We are grateful to R. Casalini, M. T. Dove, V. Heine, 
A. Navrotsky and CM. Roland for discussions, and to 
EPSRC and RFBR for support. 

APPENDIX 

In this Appendix, we establish the equivalence of d e \ 
introduced in Eq. ([T]) and the result from the viscoelas- 
tic theory of Frenkel [2(| • The relevant part of Frenkel's 
discussion starts with the modification of elasticity equa- 
tions due to relaxation process in a liquid. Consider, 
for example, the relationship between shear stress P and 
shear strain s: P = 2Gs, where G is shear modulus. In 
the presence of relaxation process the strain includes an 
extra displacement due to viscous response, and the total 
strain is written according to Maxwell interpolation as 



ds _ P 1 dP 
dt ~ + 2G~dt 

Introducing the operator 



A = l 



dt' 



(A. 14) 



(A.15) 



where r = rj/G, Eq. (|A.14|) can be written as ^§ = 
■^jAP. If A^ 1 is the reciprocal operator to A, P = 

27^- 1 §. Because f t = ^ from Eq. (|A~T5|) . P = 
2G(1 — A~ 1 )s. Comparing this with P = 2Gs, we find 
that the presence of relaxation process is equivalent to 
the substitution of G by the operator M = G(l — A^ 1 ). 

Consider the propagation of the wave of P and s with 
time dependence exp(iut). Differentiation gives multipli- 
cation by ioj. Then, ^4 = 1 + ilut, and M is: 
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M 



G 



(A.16) 



If M = Rexp(i(f)), the inverse complex velocity is 
l _ HE 



jj — v"r( cos 2 — isinj), where p is den- 
sity. P and s depend on time and position x as / = 
exp(iw(i — x/v)). Using the above expression for v, f — 
exp(iu)t) exp(— ikx) exp(— (3x), where k — LOy^cos^ 
and absorbtion coefficient f3 = w sin ^ . Combining 

27r tan ■ 



the last two expressions, we write (3 — — j — -, where 
A = =£■ is the wavelength. 

From Eq. (|A.16[) , tan <f> — -j- . For high-frequency 
waves lot ^> 1, tan<ft » = giving /? = j^. This 
is Frenkel's result [20(. Essentially the same expression 
is obtained for high-frequency longitudinal waves [20| . 
Here, Eqs. (j A. 14MA.T6]) are modified to include a finite 
zero-frequency bulk modulus. 

Lets introduce the dissipation length d = 1/(3 so that 
/ cx exp(-x/d). Then, d = Because lo = ^p, d = 

2ct (if w is close to Debye frequency, the last two formulas 
are correct approximately). Therefore, Frenkel's theory 
gives d w cr as in Eq. ([T]) . 

We note that essential to Eq. (|TJ) is the solid-like char- 
acter of LRE-induced waves: because lo 3> 1/t holds true 
for these waves, they propagate in elastic equilibrium. 
Our derivation of Eq. |T]) was based on this assertion, 
because we considered how elastic waves in a solid are 
modified in the presence of LREs. Hence, the condition 
lo ^> 1/t is used in both our and Frenkel's derivation 
above. Importantly, d c \ = cr increases with r in this 
regime of wave propagation. 

A different situation arises when lo <C 1/t. More fre- 
quently discussed in the literature, this case corresponds 
not to elastic, but to hydrodynamic equilibrium. When 
lot <C 1, Frenkel's theory gives for shear waves </> = § 
and d = Different from the high-frequency case, this 
means that low-frequency shear waves are not propagat- 
ing (because they are dissipated over the distance com- 
parable to the wavelength), a result that is also known 
from hydrodynamics [24j | . For low-frequency lot <C 1 lon- 
gitudinal waves, Frenkel's theory gives d on the order of 
— [2fl |. Here, d decreases with r, in agreement with hy- 
drodynamics 24j |. but in contrast to the high-frequency 
case due to a different regime of wave propagation. 
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